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1. Introduction 

Strange content of Nucleon ((iV|jj|iV}) and related quantities 

mAN\ss\N) dniM m, , , , , 

f T = All/ = n x _j_ ms ( N s sN ) (1.1) 

nix dm s 

has been drawing much attention in the last few years, due to its implications to the dark matter 
candidate search [1]. There have been a number of studies of (iV|jj|./V) on the lattice using different 
approaches such as direct measurement of (iV|jj|iV) matrix element with various noise reduction 
techniques for disconnected diagrams [2, 3,4], indirect measurements from fitting to ChPT formu- 
las for the nucleon mass [5, 6, 7], as well as numerical derivative via Feynman-Hellman theorem 
[8, 9]. A survey of the recent results from lattice shows there is noticeable disagreements between 
studies with relatively smaller errors while other studies with multiple lattice spacings often have ~ 
100% relative errors. This suggests the systematic errors in (JV|,k|./V) may not be well understood. 

Here we present the result of (A^|^y|A^) measurements from numerical derivative via reweight- 
ing [10, 11]. Strange quark reweighting has been successfully used to eliminate the systematic error 
from the discrepancy between dynamical strange quark mass (m s ) and the physical strange quark 
mass, and has shown to be usable in shifting m s by ~ 20% for the lattice volumes studied [12, 13]. 



2. Measurement Details 

Table 1 shows the details of the ensembles and measurements used in this study. We used (2+1) 
dynamical flavor DWF with Iwasaki gauge action generated by RBC/UKQCD Collaborations with 
a ~ 0.08, 0.1 lfm [12], which we will refer to as DWF+1 0.08fm and 0.1 lfm ensemble respectively. 
We also measured nucleon mass and (Af|jj|./V) on the DWF ensemble generated with Dislocation 
Suppressing Determinant Ratio(DSDR [14]) in addition to the Iwasaki gauge action (DWF+ID). 
Quantities necessary for continuum extrapolation of (A^|^|A^) such as lattice spacing and mass 
renormalization constant were calculated by a combined fit of DWF+I and DWF+ID ensemble 
described in [13]. 

For DWF+I 0. 1 lfm ensemble, it was observed in [12] that the combination of Coulomb gauge- 
fixed box source of size 16 and the point sink gives a good overlap with the ground state. We 
extended the measurement to use 8 different source positions per configuration, given by 

(x,y,z,t)= (0- 16,0— 16,0— 16,« x 8) for n=even, 

(12 - 28, 12 - 28, 12 - 28, n x 8) for n=odd, n = • • • 7 

Where periodicity is implied, so x=12-28 means x>=12 or x<=4. These measurements were binned 
per configuration, to get the most of the signal while avoiding the possibility of underestimating the 
autocorrelation. EigCG[15] was employed to decrease amount of computation needed for 8x12 = 
96 inversions per configuration, resulting in up to a factor of 3 decrease in the iteration number for 
the lightest mass m; = 0.005. 

For DWF+I 0.08fm ensemble, we used propagators generated by LHPC with a gaussian source 
for their study of nucleon matrix elements[16]. As it is shown in Figure 3, the nucleon propagator 
to the point sink exhibits a very slow convergence to the asymptotic value for the nucleon mass and 
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this proved to be particularly problematic when combined with strange quark mass reweighting. 
To circumvent this, propagators with unitary sink, which give consistent values and similar errors 
for unreweighted nucleon masses with point sinks but with smaller ?„,;„, were used. 

Figure 4 shows the effective mass and reweighted mass in DWF+ID (a~ l ~ 0. 14fm) ensemble. 
It is worthwhile to note that effective masses from both point and unitary sink shows a very good 
agreement, in contrast to DWF+I 0.08fm. 

Figure 1 shows the nucleon mass as a function of pseudoscalar masses squared in physical 
units. While the lack of measurements with other quantities such as the masses of other octet 
baryons does not allow us to fit the data with ChPT-motivated fitting forms with any confidence, 
the smallness of differences in masses between different ensembles suggests a 2 dependence of the 
nucleon mass is relatively small. 



ami 


am n 


MD units 


measurements 


tmin tmax 




DWF+I 0.1 lfm: V = 24 3 x 64 x 16, a" 1 = 1.75(3) Gev, am res ~ 0.003, 
am' s = 0.032, 0.033, • • • 0.04, 0.041 , • • • 0.047 


0.005 
0.010 
0.020 


0.1891(10) 
0.2420(7) 
0.322(1) 


1420,1460- •• 8980 
1460,1500- •• 8540 
1900,1920- •• 3600 


1520 
1424 
680 


5-12 


0.652(2) 
0.703(2) 
0.794(3) 


DWF+I 0.08fm: V = 32 3 x 64 x 16, a" 1 = 2.31(4)Gev 
am' s = 0.0265,0.0270- • -0.03. 


, am res ~ 0.00067, 


0.004 
0.006 
0.008 


0.1267(4) 
0.1509(3) 
0.1725(5) 


590,600- -6600 
544,552- • -7600 
590,600- •• 6600 


(1996) 
3528 
2064 


7-15 


0.474(3) 
0.503(2) 
0.524(2) 




DWF+ID: V 


= 32^ x 64 x 32, a" 1 = 1.37(l)Gev, am res ~0.001 
am' s = 0.045, 0.0455 - 0.05. 


8, 


0.01 
0.042 


0.1250(2) 
0.1810(2) 


500,508- • -2396 
608,616- •• 1920 


(904) 
1320 


5-11 


0.718(6) 
0.769(5) 



Table 1: Measurement details of ensembles used in the analysis. Parentheses in the number of measurements 
denotes the measurements are not complete. Underlines for am' s denotes the simulated strange quark mass. 



3. (A^l^^lTY) via reweighting 

Nucleon strange content is calculated by applying Feynman-Hellman theorem 

^ m = ^Ml (3 .d 

to nucleon masses for reweighted sea strange mass M^{m' s ), which is calculated by fitting the effec- 
tive masses calculated from reweighted nucleon propagator P(t,m' s ) to a constant for t = t m \ n • • • t max 
for each ensemble. 

P(t,m' s ) is calculated by taking the weighted average of the nucleon propagator P(t,m s ) with 
the weight W(m' s ,m s ), which is calculated by multiplying reweighting factors for smaller steps in 
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Figure 1: Nucleon masses of the ensembles used in this analysis in physical units. 



mass 



(P(t,m s )) = — — - r , Wi(m s ,m s ) =Wi(m s = m n ,m n -i)---Wi(m 2 ,mi =m s ). 

LiWi{m' s ,m s ) 

Using smaller steps decreases the stochastic noise in the reweighting factor, and also avoids pos- 
sible difficulties with non-gaussian distribution of the estimated reweighting factor. Reweighting 
factors Wi are calculated by estimating the ratio of strange quark determinants with Gaussian ran- 
dom vector £ . 

a i {m' n ,m n )=D[r Y D i D]{D'l)-\D i =D{[U l ],mi,m n ),D' i = {\U^m u m n ), 



Wi(m' n ,m n ) = det(Cli(m' n ,m n )y 
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(3.2) 



y/Cli(m' n ,m n ) is calculated using rational approximation, similar to what is used for the simulation 
of strange quark in [12, 13]. Using rational approximation makes Eq. (3.2) unbiased and allows us 
to use relatively noisy estimates of the reweighting factors for each mass steps. 

(Af|js|Af) for each ensemble is calculated by correlated fits to 



M^(m' s ,mi,a) = c' Q + (N\ss\N)(mi,a)m' s , 
then the continuum and chiral extrapolation is taken by fitting to 

(N\ss\N)(mi,a) = CQ + c\mi(+C2a 2 ). 



(3.3) 



(3.4) 



We used fits both with and without c^o? in Eq. (3.4), to estimate systematic error from the contin- 
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uum extrapolation. 
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Figure 2: Nucleon effective mass and reweighted mass for DWF+I 0.1 lfm ensembles 
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Figure 3: Nucleon effective mass and reweighted mass for DWF+I 0.08fm ensembles 

Figures 2 to 4 shows the reweighted nucleon mass Mjv(m^), as well as the fitting range used in 
extracting (N\ ss\N) (mi , a) for DWF+I and DWF+ID ensembles. A preliminary extrapolation to the 
continuum limit and physical pion mass gives (A^|jj|A^) = 0.33(31) at 2Gev with MS while fitting 
without a 2 gives 0.09(16), which in turn gives 



f Ts = 0.035 (33) (with a 2 ), 0.009(17) (wit out a 2 ). 



(3.5) 



4. Discussion 

We reported on the mass and the strange quark content of the nucleon, measured on the (2+1)- 
flavor dynamical DWF ensembles generated by RBC/UKQCD collaborations. Having DWF+I 
ensembles with 2 different lattice spacings allows us to do a continuum extrapolation. Comparing 
the continuum value from fits with and without a 2 terms suggests the lattice spacing error could 
be larger than previously estimated, and could explain the apparent discrepancy between different 



5 



Nucleon mass and strange content from (2+1 )-flavor Domain Wall Fermion 



Chulwoo Jung 




O O.OOl(poinl) 

- □ 0.001 (unitary) 

O 0.0042(point) 

□ 0.0042(unilary) 



' O 0.001 
- O 0.0042 



33 I 



Figure 4: Nucleon effective mass and reweighted mass for DWF+ID ensembles 
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Figure 5: Nucleon strange contents calculated via Feynman-Hellman theorem and mass reweighting, renor- 
malized at 2Gev. 



lattice studies. While we did generate new propagators on DWF+I 0.1 lfm ensembles, this results 
presented in this paper were obtained by mostly re-using propagators and reweighting factors gen- 
erated for other studies [12, 13, 16]. 

We also analyzed DWF+ID ensembles generated with Dislocation Suppressing Determinant 
Ratio(DSDR), which currently has the lightest pions and largest volumes for dynamical DWF en- 
sembles. The results with the currently available statistics are still too noisy and the results were 
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not included in the continuum extrapolation of (N\ss\N). 

It appears quite possible that further optimization of sources for the measurement of nucleon 
mass as well as increased statistics can improve the signal for (A^|^|A^) significantly without the 
need to extend the ensemble. While generating propagators with different sources and/or positions 
traditionally has required completely separate inversions and made the numerical cost prohibitively 
expensive, recent developments of various techniques such as EigCG [15]. Low/All mode averag- 
ing(LMA/AMA) [17] makes it possible to generate a large number of propagators per configu- 
rations at a numerical cost only a few times the cost for single propagator. Nucleon studies are 
particularly well suited to take advantage of this, as the correlation length for nucleons are much 
shorter than those for the light mesons. An exploratory study with AMA is under way. 

C.J was supported by the US DOE under contract DE-AC02-98CH10886. The gauge configu- 
rations used in this study were generated at QCDOC machines at Brookhaven National Laboratory 
and Edinburgh Parallel Computing Center(EPCC), and IBM BG/P machines at Argonne Leader- 
ship Computing Facility (ALCF). The reweighting factors and propagators were generated at ALCF, 
RIKEN Integrated Cluster of Clusters(RICC) and NSF Teragrid/XSEDE facilities. 
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